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Abstract. This paper presents the third in a series of single-dish studies of molecular abundances in the envelopes around a 
large sample of 18 low-mass pre- and protostellar objects. It focuses on typical grain mantle products and organic molecules, 
including H 2 CO, CH 3 OH and CH 3 CN. With a few exceptions, all H 2 CO lines can be fit by constant abundances of 7x10 **- 
8x10“^ throughout the envelopes if ortho- and para lines are considered independently. The current observational dataset does 
not require a large H 2 CO abundance enhancement in the inner warm regions, but this can also not be ruled out. Through 
comparison of the H 2 CO abundances of the entire sample, the H 2 CO ortho-para ratio is constrained to be L6±0.3 consistent 
with thermalization on grains at temperatures of 10-15 K. The H 2 CO abundances can be related to the empirical chemical 
network established on the basis of our previously reported survey of other species and is found to be closely correlated with 
that of the nitrogen-bearing molecules. These correlations reflect the freeze-out of molecules at low temperatures and high 
densities, with the constant H 2 CO abundance being a measure of the size of the freeze-out zone. An improved fit to the data is 
obtained with a ‘drop’ abundance structure in which the abundance is typically a few x lO^*** when the temperature is lower than 
the evaporation temperature and the density high enough so that the timescale for depletion is less than the lifetime of the core. 

The location of the freeze-out zone is constrained from CO observations. Outside the freeze-out zone, the H 2 CO abundance 
is typically a few xlO^'* - 10“*. The observations show that the CH 3 OH lines are significantly broader than the H 2 CO lines, 
indicating that they probe kinematically distinct regions. CH 3 OH is moreover only detected toward a handful of sources and 
CH 3 CN toward only one, NGC 1333-IRAS2. For NGC 1333-IRAS2, CH 3 OH and CH 3 CN abundance enhancements of two- 
three orders of magnitude at temperatures higher than 90 K are derived. In contrast, the NGC 1333-IRAS4A and IRAS4B 
CH 3 OH data are fitted with a constant abundance and an abundance enhancement at a lower temperature of 30 K, respectively. 

This is consistent with a scenario where CH3OH probes the action of compact outflows on the envelopes, which is further 
supported by comparison to high frequency, high excitation CS J = 10-9 and HDO line profiles which uniquely probe warm, 
dense gas. The extent to which the outflow dominates the abundance enhancements compared with the passively heated inner 
envelope depends on the filling factors of the two components in the observing beam. 
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1. Introduction 

The chemistry of organic molecules in the envelopes around 
low-mass protostars is likely to reflect directly in the molecu¬ 
lar composition of their circumstellar disks and eventual pro¬ 
toplanetary systems. A number of competing mechanisms are 
important in regulating the chemistry in these early deeply em¬ 
bedded stages: the heating of protostellar cores due to central, 
newly formed stars results in evaporation of ices in the in¬ 
nermost regions whereas shocks related to the ubiquitous out¬ 
flows may liberate ice mantles and trigger similar effects but 
on larger scales. These mechanisms have been suggested to be 
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the cause of enhancements of, e.g., H ^O, H?CO and CH 3 OH 
on small scales in the envelopes (e.g.. ICeccarelli et all Il99id 
l200(ALSch6ier et all 12002112004 . This paper, the third in ; 


a se¬ 
ries, presents an analysis of, in particular, H 2 CO and CH 3 OH 
abundances in the sample of 18 protostars studied in a wid e 
range of other molecules by LlOrgensen et alJ ilOO'A l2004d) . 
Those papers discussed observations of molecular species pre¬ 
dominantly probing the outer cold envelopes around these ob¬ 
jects. It was found that freeze-out at low temperatures and high 
densities dominates the chemistry and that, in particular, the 
freeze-out of CO is reflected in the abundances of a number 
of related species at large distances from the central protostar. 
This paper compleme nts the study of H 2 CO a nd CH 3 OH in 
a subset of objects bv lMaret~ IiT il2nn4liini . Both species 
are typical grain-mantle products observed in interstellar ices. 
To fully appreciate their chemistry it is important to compare 
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their abundances with the more general c hemical network. For 
example, high resolution observations of ISchoier et al.l (l2004l) 
indicate that the H 2 CO abundance structures may be related 
to the “drop abundance” stru ctures inferred from CO observa¬ 
tions (|j 0 rgensen et al.LEoOSi) . In these drop abundance profiles, 
freeze-out occurs in a limited region of the envelope where 
the temperature is low enough to prevent immediate desorp¬ 
tion (typically < 40 K) but the density still high enough that the 
timescale for depletion is less than the lifetime of the core (see 
also Fig. 10 of ISchdier et al.l ll2004lo . 

The class 0 protostar, IRAS 16293-2422, has long been the 
template for astrochemistry studies of dee ply embedded low- 
mass protostars du e to its rich spectrum llBlake et in 11994 
Ivan Dishoeck et all 1 19951) . IRAS 16293-2422 has a central 
warm and dense g as core where ices evaporate. Recently, 
ICazaux et alJ ll200.3l) have shown the existence of a large num¬ 
ber of complex organic species in IRAS 16293-2422, further 
underscoring the rich chemistry of this particular source. It re¬ 
mains an interesting question whether this simply reflect “first 
generation” evaporation of organic molecules at high temper¬ 
atures or whether the timescales are indeed long enough that 
a “second generation” hot core chemistry can evolve in the 
innermost regions o f these envelopes (see, e.g., discussion in 
ISchoier et al.Ll2003) . 


To put the IRAS 16293-2422 results in context, it is 
important to expand th e sample of well-studied protostars. 
|j 0 rgensen et iD il20n4cl) presented a survey of molecular 
species probing the cold outer component of protostellar ob¬ 
jects with different envelope masses, i.e., both “class 0 ” and 
“class I” objects. Although large variations in abundances oc¬ 
cur within the sample, it was found that IRAS 16293-2422 is 
in no way unique. The same objects were observed in transi¬ 
tions of H 2 CO and CH 3 OH at the JCMT iMaret et alJ J2064 
presented the H 2 CO observations, together with observations 
from the IRAM 30 m telescope, for a subset of exclusively 
class 0 objects. They reported the existence (or possibility) of 
H 2 CO abundance jumps (i.e., abundance enhancements in the 
innermost regions of the envelopes at T > 90 K where all 
ices evap orate), in some cases up to four orders of magnitude. 
However, ISchoier et alJ ll2004l) showed through high angular 
resolution data of IRAS 16293-2422 and L1448-C, that the ex¬ 
act abundance structure of the outer envelopes may severely 
affect the interpretation of the innermost envelope. For these 
low luminosity sources the warm inner regions have diameters 
<100 AU (< 0.5"), i e., are significantly diluted for single-dish 
ob servations with typical b eam sizes of 15 - 30". 

Buckle & Fullerl 1 I 2 OO 2 I) studied the low excitation (3k-'2.k) 


lines of CH3OH toward a large sample of class 0 and I objects. 
They found that a large fraction of the sources, predominantly 
the class 0 objects, show lines with two velocity components 
with CH3OH being en hanced by up to two orders of magni¬ 
tude. iBiickle & Fullerl suggested that the broad component is 
due to outflow generated shocks heating the envelope material 
and thus liberating the grain mantles. Similar effects are also 
observed in several well-studied “isolat ed” outflows well sep¬ 
arated from the protostars t h emsel v es (iBachiller et all Il995l: 
Bachiller & Perez Gutierrej, Il997t .l 0 rgensen et al.i 2004 al) . 


These and other studies illustrate that for CH3OH, a big is¬ 


sue may be whether the abundance enhancements derived from 
single-dish observations toward protostellar cores are related to 
passive heating or the action of outflows. 

This paper expands the work of iMaret et all (120041200.51) 
through observations and modeling of H2CO, CH3OH and 
CH3CN emission fo r the entire sample of pre- and protostel¬ 
lar cores studied by Ll 0 rgensen et al.l 1 I 2 OO 2 I l2004cl) . adopting 
the same physical models and approach as in these papers. In 
addition, we present high excitation CS J - 10-9 and HDO ob¬ 
servations which uniquely probe the dense and warm gas in the 
envelope. Sect.|2lpresents the observations forming the basis of 
this study. Sect. 0 describes the model approach, highlighting 
the sim ilarities and differences with the work of IMaret et alJ 
(l2004l) . Sect. 13 discusses the results, focusing on the distinc¬ 
tion between passively heated and shock processed material. 
Based on this study, it is discussed which objects are good can¬ 
didates for further studies of low-mass protostars with a “hot 
core” chemistry. 


2. Observations 

2.1. General issues 

The sample of 18 pre- and protostars presented by 
l.l 0 rgensen et al.l 1 I 2 OO 2 I) was observed in a wide range of lines at 
the James Clerk Maxwell Telescope from 2001 through 2003. 
The observations of othe r species besides H 2 CO, CH 3 OH and 
CH 3 CN are discussed bv l.l 0 rgensen et al.) (l2004cl) . iMaret et alJ 
(I20()4ll20()^ reported H 2 CO and CH 3 OH observations, respec¬ 
tively, toward the sources with the most massive envelopes - 
and hence strongest lines - in the sample, which are also in¬ 
cluded in this paper. The following sections discuss the obser¬ 
vations of each of these species in detail. All lines were ob¬ 
served and reduced in a standard way; pointing was checked 
regularly at the telescope and typically found to be accurate to 
within a few arcseconds. Calibration was checked by observa¬ 
tions of line standards and found to be accurate to within 20 %. 
The A3 and B3 receivers at 1.3 and 0.8 mm were used: the tele¬ 
scope beam sizes are typically 21" and 14" at these frequen¬ 
cies. The velocity resolution ranged from 0.13 to 0.55 km s ' 
for the different line settings. Low order polynomials were sub¬ 
tracted and the spectra were brought onto the Lmb scale through 
division by the main beam efficiencies, //mb- Values of //mb of 
0.69 and 0.63 were adopted for the 230 GHz observations with 
receiver A3 and the 345 GHz observations with receiver B3, 
respectively. 

Most of the observed lines are remarkably symmetric and 
for these lines, a single Gaussian could be fitted. The only ex¬ 
ceptions are a few of the lowest excitation lines, which are inte¬ 
grated over +2 km s“' from the systemic velocity. This velocity 
interval covers the emission from the quiescent envelope mate¬ 
rial as judged from the optically thin Gaussian lines. The line 
intensities for all species are given in Tables The corre¬ 
sponding line widths are given in Tables IA.lHA.2l in the ap¬ 
pendix. For the non-detections, 3cr upper limits are reported 
with cr - 1.2 V'JvAovT’rms where Sv is the velocity resolution, 
Aqv the expected line width to zero intensity (assumed to be 
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4 km s '), ^RMS the RMS noise level for the given resolution 
and the factor 1.2 introducing the 20 % calibration uncertainty. 

In addition to these observations, the high frequency RxW 
receiver was used to observe CS 7 = 10 - 9 at 489.751 GHz 
for four sources (L1448-C, NGC 1333-IRAS2, -IRAS4A and 
-IRAS4B) over two nights in November 2002. Special care 
was taken with the calibration: comparison with nearby spec¬ 
tral standards was found to vary by < 20 % over these two 
nights, during which the sky opacity was < 0.05 at 225 GHz 
and the elevation of the sources higher than 50°. Still, high 
system temperatures of up to ~ 5000 K were found and this, 
together with the pointing uncertainties of a few arcseconds of 
the JCMT (compared to a beam size of 10"), may cause the 
absolute calibration to be somewhat uncertain for these obser¬ 
vations. The CS 10-9 lines were detected toward NGC 1333- 
IRAS4A and IRAS4B but not L1448-C and NGC 1333-IRAS2. 
The implications of these observations are discussed further in 
Sect. 10 

2.2. H 2 CO 

In addition to the observations presented bv iMaret et alJ ll2004ll . 
H 2 CO emission from the ortho 5 15 - 4 i 4 line at 351.768 GHz 
and the para 5 o 5 - 4 o 4 line at 362.736 GHz was observed for 
all sources. Furthermore JCMT archival data exist for the para 
3 o 3 - 2 o 2 and 322 - 221 lines at 218.222 and 218.475 GHz for 
most sources and these observations wer e included for thos e 
sources not observed at the IRAM 30m bv iMaret et al.l ll2004l) . 
Table [TJlists the line intensities for all sources. The typical line 
widths of the H 2 CO lines are 1-2 km s * (FWHM). For the 
prestellar cores, LI544 and L1689B, intensiti es of 7 - 2 - 1 
and 3- 2 lines from IRAM 30 m observations bv iBacmann et alJ 
1 I 2 OO 3 I) were included in the modeling. 

2.3. CH 3 OH 

The CH 3 OH Ik -6k observations o f the sources in NGC 1333 
are presented in iMaret et al.l il200.5l) . In addition, the 5k - 4-k 
band at 241 GHz was observed for the NGC 1333 sources. 
The line intensities from the 5 a- - 4a: band data are given in 
Table El Table Ellists the CH 3 OH observations for the remain¬ 
ing sources. Only the NGC 1333 sources have detections for 
a large number of CH 3 OH lines. Besides these, only L723, 
L1448-C and VLA 1623 show detections in the Ik - 6k band 
and only of 1-2 lines each. The HDO 2ii - 2 i 2 transition at 
241.5616 GHz was covered in the CH 3 OH 5Ar - 4 a: setting, but 
was only detected toward NGC 1333-IRAS2. 

As for H 2 CO, the CH 3 OH lines can be divided into two 
forms. A- and E-type, depending on the rotation of the methyl 
CH 3 -group with respect to the OH group. Some systemat- 
ics, which may give clues to the origin of the CH 3 OH emis¬ 
sion, can be directly inferred from these results. For exam¬ 
ple, NGC 1333-IRAS2 shows the highest excitation lines for 
both the A and E-type for the 5k-4-k and Ik - 6k transitions 
compared to NGC 1333-IRAS4A and -IRAS4B. Still, the low¬ 
est excitation lines toward NGC 1333-IRAS2 are weaker than 
those toward the NGC 1333-IRAS4 sources: barring pure exci- 


Table 1. Integrated H 2 CO line intensities (f Tmb dr 
[Kkms-']). 



3o3 “ 202 * 

P-H 2 CO 
322 “ 221 * 

5o5 “ 4o4 

0 -H 2 CO 
5i5 - 4|4 

L1448-I2 



<0.9 

1.2 

L1448-C" 

3.4 

0.4 

1.3 

1.0 

N1333-I2" 

4.9 

1.0 

1.8 

1.6 

N1333-I4A" 

9.3 

2.2 

2.9 

5.5 

N1333-I4B" 

9.6 

4.7 

5.9 

7.5 

L1527" 

3.0 

0.2 

0.4 

1.0 

VLA 1623" 

5.0 


1.2 

0.9 

L483 

2.3 

< 0.2 

1.2 

1.3 

L723 

1.1" 

< 0.1" 

1.1 

2.0 

LI 157" 

1.1 

< 0.3 

0.5 

1.2 

CB244 

0.9 

< 0.4 

0.6 

1.4 

L1489 

<0.3 

<0.3 

<0.4 

0.7 

TMRl 

0.6 

<0.1 

<0.3 

0.4 

LI 544 

0.3'' 


<0.5 

<0.4 

L1689B 

1.7'' 

1.0'' 

<0.3 

<0.4 


"Observations previously reported bv lMaret~ ail <2004 . For some of 
these sources additional 0 -H 2 CO 2i2 - 1 u and 4i4 -Sn, P -H 2 CO 524 - 
423 and o-H^^CO observations and limits were reportedb^jMara^^lJ 
<2004 . *’3 o3 - 2 o 2 and 322 - 221 line intensities from lMarete^iIl2004 
are from IRAM 30 m observations. "From SEST (HPBW=24") o bser- 
vations. “^IRAM 30 m observations from iBacmann et iO <2003h : for 
these two pre-stellar cores 0 -H 2 CO 2,2 - lu observations were also 
reported hv iBacmann et alj and included in the modeling. 


tation/opacity effects, this could be interpreted as a warmer in¬ 
terior, possibly with higher abundances in NGC 1333-IRAS2. 
In particular, since the sources otherwise have relatively similar 
physical structures these differences suggest that one or more 
sources will show chemical gradients. Furthermore the CH 3 OH 
lines are significantly broader with line widths of 4-6 km s ' 
(FWHM) com pared to 1-2 krn s ' found fo r H 2 CO and the 
other species in llOrgensen et 'Zl (120021 Entnl . 

2.4. CHjCN 

CH 3 CN 14a- - 13a: was observed at 257.5 GHz for the 
NGC 1333 sources with 2 hours integration time per source, 
reaching RMS levels of 0.02 K (Tmb) in 0.36 km s ' chan¬ 
nels. At these levels the K = 0,... ,3 components are detected 
for NGC 1333-IRAS2, as shown in Fig.[I] Since NGC 1333- 
IRAS2 typically shows weaker lines for othe r molecules than 
the two NGC 1333-IRAS4 sources (see, e.g.. l.IOrgensen et ^ 
l2004cl) . this suggests that CH 3 CN may probe a different chemi¬ 
cal regime (together with the CH 3 O H lines) than the bulk, cold 
envelope in this source. Recently iBottinelli et alJ ll2004l) re¬ 
ported detections at the 40-70 mK (7’MB(peak)) of the CH 3 CN 
\Ak - 13a- components from IRAM 30 m observations toward 
NGC 1333-IRAS4A. These are consistent with our upper lim¬ 
its taking into account the smaller beam size of the IRAM 30 m 
telescope. 
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Table 2. CH 3 OH 5A:-4jf line intensities (J Tme dv [Kkm s ']) 
for sources in NGC 1333. 


Line 

Erequency 

IRAS2 

IRAS4A 

IRAS4B 


— 

4a- band; E-type 


-hOE 

241.7002 

0.41 

2.9 

2.1 

-IE 

.7672 

0.86 

7.8 

4.6 

-4E 

.8132 

0.17 

<0.06 

<0.06 

+AE 

.8296 

0.16 

<0.06 

<0.06 

-r3E 

.8430 

0.30 

0.48 

0.55 t 

-3E 

.8523 

0.17 

<0.06 

<0.06 

-rlE 

.8790 

0.33 

1.3 

1.3 

±2E 

.9044 

0.45 

2.7 

2.3 


- 

4a: band; A-type 


-hOA 

241.7914 

1.1 

8.9 

5.3 

±4A 

.8065 

0.18 

<0.06 

<0.06 

±3A 

.8329 

0.33 

0.64 

0.53 

-2A 

.8430 

0.30 

0.48 

0.55 t 

+2A 

.8877 

0.23 

0.49 

0.50 


“t”The 5-4 +3E and -2A lines are blended at 241.8430 GHz and 
have therefore (although observed) been excluded from the modeling. 
The quoted intensity refers to the total intensity of both lines. 

Table 3 . CH 3 OH line intensities (f TMedv [K km s ']) and 
limits for sources not in NGC 1333. 


7-6 


-IE OA-h 

LI 448-12 

< 0.09 

L1448-C 

0.254 0.415 

L1527 

< 0.12 

VLA1623 

0.063 0.060 

L483 

< 0.18 

L723 

< 0.06 0.373 

LI 157 

< 0.18 

CB244 

< 0.12 

L1489 

< 0.09 

TMRl 

< 0.09 



Freq. [GHz] 


Fig. 1. CH 3 CN 14a: - 13a' line observations at 257.5 GHz of 
the NGC 1333 sources. The vertical lines indicate the expected 
locations of the K = 0,1,2,3 lines. The dashed line indicates 
the 3cr detection limit. 

3. Modeling 

To model the chemica l abunda nces the approach describ ed in 
Ijdrgensen et alJ ilOO'A l2004d) and ISchoier et alJ (l2002l) . and 
utilized for the entire sample of sources and molecules, was 


Table 4 . CH3CN and CH3OCH3 line intensities (f Tmb dv 
[K km s“']) and 3cr upper limits for the sources in NGC 1333. 


Line 

Frequency 

IRAS2 

IRAS4A 

IRAS4B 

CH 3 CN 

143-133 

257.4828 

0.20 

<0.09" 

< 0.1 

142-132 

257.5076 

0.10 

- 

- 

14i-13, 

257.5224 

0.17 

- 

- 

14o-13o 

257.5274 

0.15 

- 

- 

CH 3 OCH 3 

13i,i3-12o,i2 

241.9468 

<0.08 

< 0.1 

< 0.1 


^Bottinelli et alj l2004h report detection of the 14 a:-13a- lines at the 
40-70 mK level from IRAM 30m observations toward NGC 1333- 
IRAS4A. 


adopted. Each spe cies was modele d with the envelope phys¬ 
ical structure from Lleirgensen et'al] 1 I 2 OO 2 I) derived from dust 
radiative transfer modeling of their submillimeter (SCUBA) 
continuum emission and SEDs. The line radiative transfe r was 
then performed using the code of ISchoier et al.l ll2002h con¬ 
straining the average molecular abundances. This code was 
benchmarked to high accuracy against a large number of other 
line radiative transfer cod es for a number of test problems by 
Ivan Zadelhoff et °i 1 il 2 nn 2 ll and found to agree within the Monte 
Carlo noise. 

A number of different trial abundance profiles are consid¬ 
ered besides prohles with a constant ab undance, An, of a give n 
molecule at all radii (see also Eig. 10 of ISchoier et alJ ll2004l) ). 
In “jump” prohles the abundance increase from Xq to Xj at a 
given temperature in the envelope, e.g., simulating ice man¬ 
tle evaporation. Likewise “anti-jump” prohles have a decrease 
in abundance in the inner region. Einally, in “drop” prohles the 
abundance Xq is high (comparable to the undepleted abundance 
expected in a warm molecular cloud without signihcant freeze- 
out), in the inner and outermost regions of the envelope. At 
radii where the temperature is lower than a molecule-specihc 
desorption temperature T^v and the density higher than the den¬ 
sity Wde where the timescale for freeze-out is lower than the 
lifetime of the core, the abundance decreases to Xq. 

Eor H 3 C O the collisional rate coefficients used in 
Schoier et alJ ll2003) were adopted. As described by 
Maret et al.l ll2004l) . ortho H]^CO lines are detected for 
the three sources in NGC 1333. The corresponding abun¬ 
dances were likewise calculated and taken into account in 
subsequent dis cussions. Eor CH 3 O H, new collisional rate 
coefficients by IPottage et al.l ll2004l) were used. Eor CH 3 CN 
LTE excitation was assumed. The molecular data hies are 
taken from the Lei den Atomic and Molecular Database 
dSchoier et all l2005l) and are publically available'. Each of 
the ortho- and para H 2 CO and the A- and E-type CH 3 OH are 
treated as separate molecules, which is possible since radiative 
transitions between the different species are ruled out. In the 
hrst iteration, the abundances are kept constant throughout 
the envelope and the results from the best hts are given in 
Tables ISfTl below. 

' http://WWW.strw.leidenuniv.nl/~moldata 
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Table 5. Inferred abundances for 0 -H 2 CO, P-H 2 CO and o-Hj^CO, reduced and number of lines where applicable. 


Source 

A/env 

0-H2CO 

P-H2CO 



Me 

Abundance "hms 

Abundance x^d 

^lines 




Class 0 objects 

> 0.5Mo) 




L1448-12 

1.5 

3.0x10“"’ 


1 

<1.8x10“"’ 


(1) 

L1448-C 

0.93 

8.9x10“"’ 

1.2 

3 

6.8x10“"’ 

1.0 

4 

N1333-I2 

1.7 

4.3x10“"’ 

0.63 

3 

3.9x10“"’ 

1.3 

4 

N1333-I4A 

2.3 

3.4x10“"’ 

0.59 

3 

2.3x10“"’ 

3.1 

4 

N1333-I4B 

2.0 

8.7x10“"’ 

6.2 

3 

4.1x10“*’ 

0.96 

4 

LI 527 

0.91 

8.4x10“"’ 

7.7 

3 

8.4x10“"’ 

0.44 

3 

VLA1623 

0.22 

7.7x10“"’ 

0.34 

3 

1.3x10“*’ 

0.061 

2 

L483 

4.4 

6.3x10“"’ 


1 

3.3x10“"’ 

6.7 

2 

L723 

0.62 

3.1x10“*’ 


1 

9.6x10“"’ 

2.6 

2 

L1157 

1.6 

9.2x10“" 

3.9 

3 

3.6x10“" 

2.1 

2 



Class I objects (Meny 

< 0.5Mo) 




CB244 

0.28 

5.3x10“*’ 


1 

6.9x10“"’ 

10.9 

2 

L1489 

0.097 

2.7x 10“*’ 


1 

<1.4x10“*’ 


(1) 

TMRl 

0.12 

5.0x10“*’ 


1 

3.4x10“*’ 


1+(1) 



Pre-stellar cores 




LI544“ 

2.5 

4.0x10“" 


1 

3.0x10“" 


1 

L1689B“ 

2.5 

2.0x10“"’ 

8.4 

2+(l) 

1.2x10“"’ 

0.6 

2 

0 

n 

0 

N1333-I2 

1.7 

3.0x10“" 

3.3 

2 




N1333-I4A 

2.3 

5.2x10“'^ 

2.5 

3 




N1333-I4B 

2.0 

8.0x10“" 


1 





“Quoted value based on line intensities reported by iBacmann et all l2003h . Upper limits from JCMT lines only are one to two orders of 
magnitude larger. 


Table 7. Inferred abundances and reduced from models of 
CH 3 CN toward IRAS2. 


Source 

Abundance 


^lines 

N1333-I2 

CH3CN 

8.3x10“" 

8.5 

4 


7x10“*’“ 

1.4 

4 


“Abundance in inner {J > 90 K) region. 3cr upper limit on the abun¬ 
dance in the outer envelope of 2x 10“'*. 


3.1. H 2 CO 

For most sources the H 2 CO line intensities are well-fit with 
constant abundance models for each of the P-H 2 CO and o- 
H 2 CO species. For a few sources, one of these two has a 
higher than 3 ( 0 -H 2 CO; NGC 1333-IRAS4B, L1527, 
lTi 57 and L1689B; P-H 7 CO : NGC 1333 - IRAS 4A, L483 and 
CB244). This is in contrast to lMaret et alJ ( l2004t) who infetTed 
large abundance jumps for the studied sources. As shown be¬ 
low this is due to a number of d i fferen ces in the assumptions 
between the work of iMaret et alJ ( l2004l) and this study, in par¬ 
ticular the ortho-para ratio. Below we discuss some of these 
differences. 


3 . 1 . 1 . Velocity structure 

|m aret et a i 1 1I2004 assume a non-turbulent but infalling enve¬ 
lope, whereas our work assumes a constant turbulent broad¬ 
ening of 0.5-1 km s ' throughout the envelopes reproducing 
the observed line widths. The derived abundances do not de¬ 
pend on the velocity held if onl y optically thin lines are con¬ 
sidered (l.lOrgensen et alll2004cll . For the constant abundances 
listed above, the observed lines have typical optical depths of 
0.1-1. This constant width may not be an adequate description 
of emission originating from the innermost dense regions in the 
case of large abundance jumps, but in the context of an inside- 
out collapsing envelope, the outer envelope will naturally be 
characterized by the turbulent broadening. In class 0 objects 
typical infetTed ages are ~ lO"* years, which for typical values 
of the sound speed of 0.3-0.5 km s“* translate to infall radii 
of 500-1000 AU. Such sizes are unresolved by single-dish ob¬ 
servations and even by medium-high resolution interferometer 
data. This implies that the infall radius encompasses the hot in¬ 
ner region, but that most of the mass in the envelope material is 
not infalling. 

From line radiative transfer it is possible to make exact pre¬ 
dictions for the line prohles, which can be compared to the ob¬ 
served spectra to constrain the velocity held. Fig. |2] compares 
the observed 5 o 5 - 4 o 4 and 5 i 5 - 4 i 4 spectra toward L1448-C with 
the modeled line prohles adopting the non-infalling but turbu¬ 
lent envelope model (with constant H 2 CO abundance) from this 
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Table 6. Inferred abundances for CH 3 OH, reduced and number of lines where applicable. 


Source 

CH 3 OH A-type 
Abundance Xka 

^^lines 

CH 3 OH E-type 
Abundance Xred 

^lines 


Class 0 objects (Mc„v 

> 0.5Mo) 



L1448-12 

<1.1x10-“* 


(1) 

<1.4x10-“* 


(1) 

L1448-C 

1.3x10-**“ 


1 

1 .0x10-“*“ 


1 

N1333-12 

1.4x10-**“ 

20.4 

8 

1.3x10-“*“ 

18.4 

13 

N1333-I4A 

2.9x10-**“ 

4.5 

7 

2.5x10-“*“ 

3.1 

7 

N1333-I4B 

3.5x10-**“ 

11.4 

7 

9.1x10-“*“ 

9.6 

7 

LI 527 

<5.5x10-“* 


(1) 

<6.8x10-*'* 


(1) 

VLA1623 

2.5x10-“* 


1 

3.4x10-“* 


1 

L483 

<2.2x10-“* 


(1) 

<2.7x10-“* 


(1) 

L723 

1 .8x10-** 


1 

<3.6x10-“* 


(1) 

L1157 

<2.7x10-“*“ 


(1) 

<3.4x10-*'*“ 


(1) 


Class I objects (Mem 

< 0.5Mo) 



CB244 

<1.0x10-** 


(1) 

<1.3x10-“* 


(1) 

L1489 

<9.9x10-“* 


(1) 

<1.2x10-“* 


(1) 

TMRl 

<2.2x10-** 


(1) 

<2.8x10-“* 


(1) 


“See also lMaret et all j2Q05l) . 


paper and t he non-turbu l ent. in falling model (with abundance 
jump) from iMaret etl^ 1I2OO4I) . Considering only the higher 
excitation H 2 CO 7 = 5-4 lines should limit confusion from the 
surrounding material. It is seen that the “turbulent” envelope 
model clearly provides the closest match to the observed line 
prohle, whereas the non-turbulent, infalling model provides too 
much broadening. It is interesting to note that the same turbu¬ 
lent broadening which is used to fit lo w-excitation CO i s otopic 
lines from the cold, outer envelope (ilOrgensen et all l 2002 i) 
works well in also explaining the observed H 2 CO line widths. 
In the cont ext of, e.g., an inside-out collapsing envelope as in 
this suggests that the infalling region of the enve¬ 
lope is small and that even the relatively high excitation H 2 CO 
lines studied in this paper are still predominantly sensitive to 
the outer envelope. It is worth re-emphasizing that both tur¬ 
bulent and non-turbulent/infalling models give a good constant 
abundance fit when the ortho and para lines are treated indepen¬ 
dently and that abundance jumps are in general not required for 
the studied sources, other than in the context of a “drop model” 
(see Sect. 1113. Naturally a combination of turbulent broad¬ 
ening and systematic motions (e.g., infall) most likely applies 
to the studied envelopes and a fully self-consistent model ex¬ 
plaining also emission of optically thick lines would have to 
take these into account. Observations of optically thick transi¬ 
tions of, e.g., HCO^ could be used to constrain such models, 
but this is beyond the scope of the present study. 



2468 2468 


V 

LSR 


[km s-i] 


Fig. 2. Spectra and modeled line profiles of the 5 o 5 - 4 o 4 (para) 
and 5 15 - 4i4 (ortho) H 2 CO lines (upper and lower panels, re¬ 
spectively) toward L1448-C in a pure turbulent, constant abun¬ 
dance envelope (this paper; left) and in a non-turbulent enve¬ 
lope wi th an abundance jum p infalling toward a 0.5 Mq central 
object dMaret et alJ (l20()4l) : right). The spectra have all been 
normalized by division with the total integrated line intensity 
to bring out more clearly the comparison between the actual 
line shapes. 


3 . 1 . 2 . H2CO ortho-para ratio 

A second difference with iMaret et al.l(l2004l) is the assumption 
of a fixed ortho/para ratio. Since para and ortho H 2 CO can be 
considered to be separate molecules, their abundances can be 
determined independently. Fig. [^compares the abundances for 
the two species. A very close coiTelation exists (Pearson cor¬ 


relation coefficient of 0.9; see also Sect. |3T3, which can be 
htted by an ortho-para ratio of 1.6+0.3. The very tight corre¬ 
lation indicates that both species probe the same region of the 
envelope and that their abundance ratios are established under 
similar conditions in all sources. 

iMaret et alJll2004l) assumed an ortho-para ratio of 3 to com¬ 
bine 0 -H 2 CO and P-H 2 CO line observations to constrain the 
abundance structure. Whereas this gives in principle fewer free 
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Fig. 3. Comparison between the abundances of the P-H 2 CO and 
0 -H 2 CO species. The class 0 objects are indicated by hlled 
diamonds, the class I objects by open diamonds and the pre- 
stellar cores by hlled squares. The class 0 objects VLA1623 
and IRAS 16293-2422 have been singled out by an open square 
and star, respectively. Sources with poor hts ^ 5) to either 
the P-H 2 CO or 0 -H 2 CO species are indicated by black crosses. 
Of these N1333-I4B has further been singled-out with the open 
circle. The statistical uncertainties from the hts including cali¬ 
bration uncertainties are 10-20%, i.e., the ertor bars are of sim¬ 
ilar size as the plot symbols. The solid line indicates the best ht 
linear correlation between the two sets of abundances (exclud¬ 
ing the poorly ht sources), corresponding to an ortho-para ratio 
of 1.6:1. The dashed line indicates the relation for an ortho-para 
ratio of 3:1. 

parameters in the modeling, one should be careful when inter¬ 
preting the results. This is clearly illustrated in Fig. |4] where 
constraints on the total H 2 CO abundance from ortho and para 
lines for L1448-C are shown. It can be seen that the overlap¬ 
ping conhdence levels depend critically on the adopted ortho¬ 
para ratio, with a high ortho-para ratio driving an increased 
abundance jump in the inner envelope up to 4 orders of mag¬ 
nitude. For both species a constant abundance model provides 
a good ht to the observed lines and the combination of the two 
suggests an ortho-para ratio of 1.6, in agreement with the con¬ 
clusion above. Fig.^compares the modeled and observed line 
intensities for L1448-C as functions of energy level of the ob¬ 
served H 2 CO transitions. No systematics are seen (for example, 
transitions from higher levels are not systematically underesti¬ 
mated in the models) and almost all observed intensities are 
well reproduced within the observational uncertainty. In gen¬ 
eral, it is difficult to constrain the ortho-para ratio for a specihc 
source because of the intrinsic uncertainty in the observations 
and modeling, including va rying ortho-pa r a ratio s with position 
and varying optical depth dSchoier et all I2OO2I) . The strength 
of our analysis is that the ratio is based on a large sample of 
sources, statistically reducing some of these uncertainties. 

The ortho-para ratio itself contains interes ting informa¬ 
tion ab out the H 2 CO formation, as discussed bv iKahane et al.l 
lll984l) . At high temperatures the ortho-para ratio approaches 
the relative statistical weights of 3:1, but thermalization at 


Fig. 4. Constraints on the H 2 CO abundances in the inner (T > 
90 K) and outer {T < 90 K) envelope (Xj and Xo, respec¬ 
tively) and effect of adopted ortho-para ratio for L1448-C. 
A non-turbulent, f ree-falling envelope has been adopted as in 
iMaret et^ ll2004l) . The solid line contours indicate the 2 and 
4cr confidence levels for ortho-H 2 CO, whereas the dashed line 
contours indicate the corresponding conhdence levels for p- 
H 2 CO. The grey scale contours indicate similar conhdence lev¬ 
els for the H 2 CO abundance combining the constraints from the 
two datasets and assuming the ortho-para ratios of 1.0, 1.6, 2.0 
and 3.0, respectively, as indicated in the top of each panel. 



Fig. 5. Comparison between observed and modeled line inten¬ 
sities as function of the energy of the upper level for each H 2 CO 
transition toward L1448-C. The observed intensities of the para 
lines are indicated by the circles, the intensities of ortho lines 
by the squares and the modeled intensities by the crosses. 

lower temperatures {T ~ 10 - 15 K) makes the two abun¬ 
dances closer to equal. As argued bv iKahane et al.t gas-phase 
thermalization is improbable: the life-times of the established 
ortho-para species are signihcantly longer than, e.g., the time it 
takes to destroy a given H?CO molecule through reactions with 
molecular ions. IKahane et al.l furthermore hnd that chemical 
reactions lead to ortho-para ratios of 3-5. In contrast, H 2 CO 
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Table 8. Pearson correlation coefficients betw een abundances 
found in this paper and abundances from IjOreensen et alJ 

M). 



P-H 2 CO 

0 -H 2 CO 

CO 

0.63 

0.70 

HCO+ 

0.46 

0.76 

CS 

0.42 

0.60 

SO 

0.13 

0.05 

HCN 

0.72 

0.63 

HNC 

0.57 

0.74 

CN 

0.57 

0.81 

HC 3 N 

0.44 

0.70 

P-H 2 CO 

0 -H 2 CO 

0.92 

0.92 


10 - 


o 

“cm 10-9 


10 - 


o 

□ 


0.1 


1.0 


Mass / M, 


can be thermalized on dust grains and released subsequently. 
For grain temperatures of 10-15 K, an ortho-para ratio of 1.5 
is in good agreement with the results from this paper. This is 
also what is found in typical dark clouds such as TMCl and 
L134N, whereas warmer regions such as the Ori on clouds show 
ortho - para ratios closer to the statis tical 3:1 ratio dKahane et al.l 
Il984l iMangum & Woottenlll993ll . The derived ortho-para ra¬ 
tio of 1.6+0.3 for the sources in this sample therefore suggests 
that the observed H 2 CO lines are predominantly sensitive to 
the outer cold envelope and that relative abundances of the or¬ 
tho and para species are established there. 


o 

u 




10 - 


10 - 


10-19 



10-9 lQ -4 

[CO] 


3.1.3. Comparison to other molecules and implications 
for abundance structures 

To quantify the relations between the abundances of the ob¬ 
served molecular species, [lorgensen et all (l2004cl) calculated 
Pearson correlation coefficients for each pair of abundances. 
The Pearson correlation coefficient, P, is a measure of how well 
a (x,y) data set is fit by a linear correlation compared to the 
spread of (x,y) points. Values of +1 indicate good correlations 
(with positive or negative slopes) whereas a value of 0 indicates 
no correlation. In our studies of other molecules, strong corre¬ 
lations (|P| > 0.7) were found between molecules for which 
relations were expected based on chemical considerations, for 
example between CO and HCO^ or between the sulfur-bearing 
species. To extend this discussion, correlation coefficients were 
calculated between the abundance s found in this paper and 
those from l.TOrgensen et ^ il2004ci) (see Table|^. 

An important conclusion is the apparent anti-correlation 
between the H 2 CO abundances and envelope mass within the 
radiu s where the temperature drops to 10 K (I.TOrgensen et al.[ 
[ 2 ^ as shown in Fig. |6l This is also reflected in the 
correlations in Table 0 between H 2 CO and molecules such 
as CO whose abundances decline wi t h increasing mass 
jjdrgensen et all l2002l) . IjOrgensen et alJ ( l2004cl l2005h sug¬ 
gested that this correlation is well explained with a drop- 
abundance profile, with the different H 2 CO abundances reflect¬ 
ing the size of the region over which freeze-out occurs. 


Fig. 6. Total H 2 CO abundance vs. envelope mass (upper panel) 
and CO abundance (lower panel). The envelope mass refers 
to the mass within the radius where the temperature drops to 
10 K as estimated from the dust radiative transfer models (see 
LlOrgensen et al.H2002l) . For objects where only the abundance 
of one of the two H 2 CO species has been constrained, an ortho¬ 
para ratio of 1.6 is assumed. As in Fig. 13 the class 0 objects are 
indicated by the class I objects by “0” and the pre-stellar 
cores by The class 0 objects VLA1623 and IRAS 16293- 
2422 have been singled out by and respectively. 


As an illustration, the drop abundance structure is tested 
for the H 2 CO lines toward NGC 1333-IRAS4A. For the H 2 CO 
lines we first constrain the depletion density, Wde. and des¬ 
orption tempe rature, Tpy, from obser vations of the CO lines 
presented by LlOrgensen et al.l (I2flfl2h . The CO data toward 
NGC 1333-IRAS4A are well-fitted with depletion by a factor 
of 50 in the region of the envelope where the density is higher 
than 6x10^ cm“^ and temperature lower than 40 K. These con¬ 
straints are used as input for the H 2 CO chemical structure, so 
that only the overall and depleted abundances, (Vq and Xu, re¬ 
spectively) are left as free parameters. The results of these fits 
are shown in Fig. 0 both ortho and para lines are consistent 
with an abundance drop of approximately an order of mag¬ 
nitude. The confidence regions for the two H 2 CO species 
agree at the Icr level assuming an ortho-para ratio of 1.6. Also 
the o-HPCO observations agree with those of the main iso- 
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Fig. 7. Constraints on H 2 CO abundances in a “drop” model for 
NGC 1333-IRAS4A, i.e., a model where the molecule is de¬ 
pleted to Xb in the regions where the density is higher than 
tide and temperature lower than constrained from CO ob¬ 
servations. The solid and dashed contours indicate the Icr, 2cr 
and 4cr confidence levels for P-H 2 CO and 0 -H 2 CO, respec¬ 
tively whereas the grey-scale contours indicate the confidence 
region for the two species combined with an ortho-para ratio 
of 1.6. The black/white dotted line indicates the best fit rela¬ 
tion from the o-Hj^CO lines. The best fit model (2fo=3xlO'®, 
2fD=4xlO“'°) combining all lines has indi¬ 

cated with the 


topic lines assuming a ^^C:^^C ratio of 70. The best fit drop 
model has an undepleted abundance Xq = 3 x 10“^ and an 
abundance in the depletion region of Xd = 4 x 10 The re¬ 
duced for this model is 1.1 for 10 fitted lines (including all 
ortho, para and Hj^CO lines). This suggests that the variations 
in H 2 CO abundances reflect, to first order, the variations due 
to freeze-out, with the chemical network subsequently regu¬ 
lating the abundances. High-resolution interferometer observa- 
tions confirm this stru cture for IRAS 16293-2422 and L1448-C 
dSchoier et alll2003) . 

Table compares the drop abundance profiles for 
NGC 1333-1RAS4A with t hose of L1448-C and IRAS 16293- 
2422 dSchoier et ali |2004 . Similar values for Xq are found 
within a factor 3 and with the H 2 CO abundance decreased by 
an order of magnitude in the freeze-out zone. Also the maxi¬ 
mum constant abundance for the entire sample is in agreement 
with this value for Xq, supporting the suggestion that the spread 
in constant H 2 CO abundances reflects the size of the freeze-out 
zone. 

A varying abundance structure is also preferred for 
NGC 1333-IRAS4B: as seen from Table|5landFig.|3 the ortho¬ 
para ratio for this particular source is less than 1. This prob¬ 
lem is not alleviated by the introduction of a “drop” profile, 
which still gives an ortho-para ratio below unity and a poor 
fit to, in particular, the ortho-H 2 CO lines. An abundance in¬ 
crease at low temperatures, however, does a better job: Fig.|3 
shows models for NGC 1333-IRAS4B with abundance jumps 


Table 9. Summary of models with varying H 2 CO abundance 
structure. 



JfD 

^0 

NGC 1333-IRAS4A 

4x10^“’ 

3x10-" 

Liqqs-C" 

1x10-" 

1x10-* 

IRAS 16293-2422“ 

3x10-*" 

1x10-* 

Sample'’ 

7x10-“ 

8x10-" 

NGC 1333-IRAS4B' 

<1x10-" 

1x10-* 


°From l^hoier et ^j2004) fitting both single-dish and interferometer 
data. ^Minimum and maximum abundances for the entire sample in 
this paper. ''Model with abundance jump at 20 K. Ad for this source 
refer to the abundance where T < 20 K and Xq to the abundance where 
r > 20 K. 


at differing temperatures. An abundance jump at T^y < 30 K 
from ~ 10 “'® to ~ 10 “* makes it possible to fit the lines with 
an ortho-para ratio above unity, and to bring the abundance in¬ 
ferred from the o-Hj^CO lines in agreement with that of the 
o-Hj^CO lines. A jump at low temperatures also significantly 
improves the best fit for the two speci es separately. This is in¬ 
teresting compared with the results of iMaret et alJ ll2004 who 
inferTed an abundance enhancement close to 4 orders of mag¬ 
nitude in NGC 1333-IRAS4B, but with a rather low quality of 
the fit (x^^d ~ suggests that the model with a jump at 

temperatures of 90-100 K is not adequate to describe the abun¬ 
dance structure for NGC 1333-IRAS4B but that other mecha¬ 
nisms such as the action of the protostellar outflow regulate the 
H 2 CO abundance for this source. 

To summarize these discussions: by fitting the ortho and 
para H 2 CO lines independently for the larger sample of 
sources, the ortho-para ratio can be constrained statistically to 
be 1.6 + 0.3. The observed cotTelations with other species from 
the survey of LlOrgensen et al.l (l2004cl) suggest that the H 2 CO 
abundances are related to the overall chemical network - pri¬ 
marily reflecting freeze-out of CO at low temperatures and high 
densities. This is further illustrated by fits to the NGC 1333- 
IRAS4A H 2 CO data, which are improved with a “drop abun¬ 
dance” profile. In these models, the H 2 CO abundance is de¬ 
creased by an order of magnitude in the cold freeze-out zone. 
No jump in abundance in the innermost {T > 90 K) region is 
needed, although NGC 1333-IRAS4B is best fit with an abun¬ 
dance increase where T > 20 - 30 K. Observations of higher 
excitation lines are needed to constrain jumps in the hot core 
region. The NGC 1333-IRAS4B data suggest that for specific 
sources, other effects, such as the impact of the outflow, may 
play a role in defining the H 2 CO abundance structures. 

3.2. CH 3 OH and CH^CN 

The upper limits for the constant CH 3 OH abundances found 
for the sources in sample range from a fewxl 0 “'° to ~ 10 “®. 
The upper limits are typically a factor of a few below the 
abundances determined from lower excitation CH 3 OH lines by 
iBiickle & Fuller! (I2OO2I) . However, since their abundances were 
calculated relative to CO, for which abundances are lower by 
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Fig. 8 . ;^^^-confidence contour plots for P-H 2 CO and 0 -H 2 CO 
abundances (upper and lower panels, respectively) toward 
NGC 1333-IRAS4B. The minimum for the given jump 
temperature is indicated in the top of each panel. The dashed 
lines indicate the constraints on the abundances from the o- 
H^^CO lines. In the Tev = 20 K and Tev = 30 K panels the best 
ht models combining the constraints for all lines (with mini¬ 
mum = 1.7 and 3.8, respectively) have been indicated by 
For the other values of Tev, the minimum> 5. 


an order of magnitude than the canonical value in the class 
0 objects, the upper lim its from this paper and the results 
of iBuckle & Fullen 1I2OO2I) are still consistent. iBuckle & Fuller] 
also found a slight decrease in CH 3 OH abundance with bolo- 
metric temperature (i.e., lower CH 3 OH abundance in the class I 
objects). This may again be a result of the abundance calculated 
relati ve to CO, since the CO abundance varies with envelope 
mass llj 0 rgensen et alll2002h. A similar t rend w as seen for the 
sulfur-bearing species jBuckle & FuB^ l2003[ see discussion 
in Jprgensen et al. 2004d). _ 

As noted by iMaret et"^ ll2005l) . the CH 3 OH data for 
NGC 1333-IRAS2 and -IRAS4B (see also Table cannot 
be modeled with constant abundances. In the best cases such 
models gw/e x\ed ~ 10 “ 20. NGC 1333-IRAS4A also shows 
mediocre hts for each of the A and E-type species with con¬ 
stant abundances, but still better than the two other sources 
(tL-3-4). 

Our analysis uses the ne w collisional rate coe fficients for 
CH 3 OH recently published bv IPottage et al.l(l2004l) . Compared 
to the old rate coefficients, the derived line intensities vary in 
certain cases by up to 50%. However, no systematic trends are 


Fig. 9. ;^f^-conhdence contours for jump models with varying 
Tev for CH 3 OH lines toward NGC 1333-IRAS2. The minimum 
Tred given in the upper right corner of each panel. 

90 K 50 K 30 K 

= 4.7 = 4.5 x\„.. = 4.8 



Oil 


10-9 


10-9 lo-ti 10-9 10-^^ 10-9 lO-'^ 

X 

J 

Fig. 10. As in Fig. |9] for CH 3 OH lines toward NGC 1333- 
IRAS4A. 


seen and therefore the derived abundance structures are un¬ 
changed. Still, this example illustrates that the derived abun¬ 
dances - especially when based on constraints from only a few 
lines - may be uncertain by up to a factor of 2 due to uncertain - 
ties in the collisional data alone (see also lSchoier et alll2005t) . 
The poor hts can be improved by includi ng evaporation of 

S ee mantles a t temperatures > 90 K llCeccarelli et all 
Schoier et all 1200^ iMaret et all 1200411200.51) . To test 
this, a step function for the abundance was introduced, with 
a jump in abundance from Xq in the exterior to Xj in the in¬ 
terior at a radius corresponding to a specihe temperature Tev. 
Models were run for - 30, 50 and 90 K for each of the 
three NGC 1333 sources. Table fTDI gives the best ht models 
and Fig. 1911 ll show the derivedconhdence plots for each of 
the temperatures and for each of the sources. They clearly show 
different behavior: NGC 1333-IRAS2 is nicely ht with a jump 
at 90 K, whereas NGC 1333-IRAS4B is much better ht with 
a jump at 30 K. For NGC 1333-IRAS4A the models suggest a 
best ht for a constant or “anti-jump” abundance structure. 

It is also noteworthy that the best hts of both A- and E- 
type species are obtained for similar abundances at the best 
htting evaporation temperatures. As for the H 2 CO ortho-para 
species, no transitions between the A- and E-type CH 3 OH lev¬ 
els are expected. The E/A-type abundance ratio, however, only 
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Table 10. Derived CH 3 OH and CH 3 CN abundances assuming abundance jumps in the inner T > regions. 



Tev 

Species 

Aj 

Ao 

aL 



CH3OH 




IRAS2 

90 K 

A-type (8) 

9x10 " 

1x10-“' 

6.6 



E-type (13) 

3x10-^ 

7x10-'“ 

2.5 

IRAS4A 

50 K 

A-type (7) 

<5x10 " 

4x10-“ 

4.5“ 



E-type (7) 

<2x10 " 

3x10-“ 

4.4“ 

IRAS4B 

30 K 

A-type (7) 

1x10^’ 

3x10-“ 

2.8 



E-type (7) 

9x10 " 

3x10-“ 

2.7 

IRAS 16293-2422* 

90 K 

A-l-E-type (23) 

1x10^’ 

6x10-“ 

1.2 



CH3CN 




IRAS2 

90 K 

A-type 

7x10^'' 

< 3x10" 

1.4 


“No strong constraints exist on the evaporation temperature for NGC 1333-IRAS4A in accordance with the conclusion that this source is well- 
fitted by a constant abundance throughout the envelope (see Fig. llOl . ^Results for IRAS 16293-2422 from ISchoier et ^ <2002h assuming the 
abundances of the A- and E-type CH3OH to be identical. 


90 K 50 K 30 K 

. = H-8 x^ . = 3.4 : X^ . = 2.8 



10-9 


10-9 10-^ 10-9 10-9 10-'^ 10-9 10-8 10-’ 10-9 

X 

j 

Fig. 11. As in Fig. |3 for CH 3 OH lines toward NGC 1333- 
IRAS4B. 


varies from 0.69, corresponding to thermalization at 10 K, to 
unity in the high temperature limit (iFriberg et all Il988ll . The 
abundances are therefore expected to be similar for the best fit 
models, which appears to be the case as illustrated in Fig. EEU 

For NGC 1333-IRAS2 similar jump models were run for 
CH 3 CN, and the best fit abundance is shown in Fig. [21 
Interestingly the CH 3 CN lines also give an abundance jump 
of approximately two orders of magnitude at 90 K, similar 
to what is found for the CH 3 OH lines. Again jumps at lower 
temperatures are not favored for this source. The value of 
X] for NGC 1333-IR AS2 is comp^able t o th at derived for 
IRAS 16293-2422 by LSchoier et alJ (l 2002 h and ICazaux et al.l 
il2003ll . There is no strong evidence to suggest that the chem¬ 
istry of this molecule is significantly different in the innermost 
envelopes around the objects discussed in this paper compared 
to IRAS 16293-2422. 


90 K 50 K 30 K 



Fig. 12. ^^-confidence contours for jump models with varying 
Tev for CH 3 CN lines toward NGC 1333-IRAS2. The minimum 
Tred given in the upper right corner of each panel. 


4. Discussion 

4.1. Hot core vs. outflow 

As discussed in Sect. im the H 2 CO abundances are found 
to be consistent with constant abundances throughout the en¬ 
velopes for most sources. The CH 3 OH results for the three 
NGC 1333 sources, in contrast, imply abundance variations, 
but these occur over significantly different regions of the en¬ 
velope. The three objects have rather similar density and tem¬ 
perature profiles and the observed differences therefore sug¬ 
gest other causes for the abundance enhancements than passive 
heating of the envelope material. 

An important clue comes from the differing line profiles of 
H 2 CO and CH 3 OH. Fig.ll3lcomDares the profiles for the H 2 CO 
5 o 5 - 4 o 4 and CH 3 OH 7_i - 6 _i-E lines toward NGC 1333- 
IRAS2. The CH 3 OH line is significantly broader, with a width 
of 4 km s ' compared to the 1.5 km s ' for the H 2 CO line. 
Part of this could be due to differences in the thermal broaden¬ 
ing if CH 3 OH probes warmer gas much deeper in the envelope. 
The radiative transfer models, however, take this explicitly into 
account and it is concluded that a significantly higher turbu¬ 
lent broadening is required to model the observed CH 3 OH lines 
compared with those of H 2 CO and other species. This suggests 
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Fig. 13. Observed and modeled spectra of H 2 CO Sqs - 4 o 4 and 
CH 3 OH 7_i6_i-E lines toward NGC 1333-IRAS2. The two 
models shown use turbulent line broadening of 0.8 km s“', 
which can also accoun t for, e.g., the C '*0 lines modeled in 
Ijdrgensen et alJ ( l 2002 t) (solid line), and 2.5 km s * (dashed 
line), respectively. Note how the two lines probe significantly 
different velocity fields in the envelope. 


that the CH 3 OH lines probe a differ ent part of the envelope th an 
H 2 CO and the species discussed bv l.Te>rgensen et alJ ll2004cl) . 

One explanation for the enhancements and profiles of 
CH 3 OH could be the impact of the outflows on the inner en¬ 
velopes. This was sugges ted to be the case for NGC 1333- 
IRAS4A and IRAS4B by iBlake et aP lll995h . who noted that 


the abundance enhancements can occur in the outflow shear 
zones probed by the very broad, kinematically distinct CS and 
CH 3 OH lines. Large enhancements of CH 3 OH are seen in the 
shocks driven by the outflows from a number of protostellar 
sources where the shock is well separated from the central 
prot ostar including L1448-C, LI 157 and NGC 1333-IRAS2 
(e.g., BachilleretaL^ 199-5l:lBachiller & Perez Gutierre 7 lll 997 t 


Jorgensen et all l20()4ai) . CH 3 OH is one of the molecules that 

shows the largest jumps in abundances between cold and warm 
gas and therefore traces more clearly the origin of the abun¬ 
dance enhancements. 

An outflow scenario could explain the relatively low des¬ 
orption temperature found for the CH 3 OH lines as well as the 
fits to the H 2 CO lines for NGC 1333-1RAS4B discussed in 
Sect. 13.1.31 If these species are desorbed from grains due to 
the action of the outflow, they could thermalize at temperatures 
closer to that of the envelope, i.e., lower than in the small “hot 
core” region. NGC 1333-1RAS4A is remarkable since the fits 
to the CH 3 OH lines do not require a jump in abundance. This 
difference from NGC 1333-1RAS4B coul d be caused by the 
differe nces in the outflow morphologies: IPi Francesco et all 
1I2OOII) imaged the NGC 1333-1RAS4 region at high resolu¬ 
tion and found that the outflows are probed by the wings of 
H 2 CO and CS. The images show that the 1RAS4B outflow 
is more compact with emission over a ^ 15" region com¬ 
parable to the 350 GHz JCMT single-dish beam. The mor¬ 
phology thus suggests that enhancements for the NGC 1333- 
1RAS4B outflow occur on scales where the single-dish obser¬ 
vations are the most sensitive. Also the SiO J =5-4 line at 
217.1 GHz is only detected toward NGC 1333-1RAS4A and 
not 1RAS4B. This indicates that the CH 3 OH enhancements in 
1RAS4B are not directly associated with a high velocity shock 
but more likely result from the shear between the envelope 




10' 10^ 10® 10^ 10' 10® 10® 10^ 

R [AU] R [AU] 


Fig. 14. Comparison of the inferred temperature, density 
and abundance structures in NGC 1333-1RAS2 (left) and 
NGC 1333-1RAS4B (right). The upper panels show the tem¬ 
perature (black solid line) and density (grey dashed line) pro¬ 
files for each source. The lower panels show the inferred CO 
drop abundance (black soli d line) and CH 3 OH jum p abundance 
(grey dashed line) profiles iltlrgensen et al.[E005[ this paper). 


and out flow. Note also that direct images of outflow induced 
shocks jBachiller & Perez Gutierre jll997HGarav et alll200(A 

lBachillereraLn200ll: J 0 reen^n et al.L 2004al) indicate that 

CH 3 OH does not survive at the highest shock speeds, which 
could further explain the differences between 1RAS4A and 
1RAS4B. In contrast, the impact of the outflow in NGC 1333- 
IRAS2 is on much smaller scales, more heavily diluted in the 
single-dish beam. Fig. m compares the inferred abundance 
structures for NGC 1333-IRAS2 and -IRAS4B. 


4 . 2 . CS J = 10-9 and HDO as dense gas probes 


Further support for the importance of the outflows comes from 
observations of high excitation CS J - 10-9 and HDO lines 
obtained with the JCMT. For CS J - 10-9, broad lines 
(FWHM ss 8 km s ') are detected toward both NGC 1333- 
1RAS4A and 1RAS4B (Fig.[T5}, which lack the central narrow 
peak seen for the lower excitation lines (Fig. d. Although, 
the absolute calibration may be somewhat uncertain, the CS 
J - 10-9 line is approximately 5 times stronger toward 
1RAS4B than 1RAS4A, consistent with the above conclusion 
that the dense outflow gas fills a larger fraction of the beam for 
1RAS4B. This is in contrast to the lo wer excitation CS lines re - 
ported, e.g., bv iBlake et al.l(ll99^ and ljprgensen et an(l2004d) . 
which support a more compact origin of the CS outflow emis¬ 
sion in 1RAS4B than 1RAS4A. 

Table rm compares the predictions for the CS J - 10-9 
line intensity assuming C S abundances in the qui escent enve¬ 
lope from the models of l.lc>rgensen et al.l ll2004di . The enve¬ 
lope models predict significantly less emission than observed 
for 1RAS4A and 1RAS4B, but the non-detection toward 1RAS2 
is con sistent within the 3 cr nois e level. Note that the CS intensi¬ 
ties in .lOrgensen et ^( l2004cl) were found by integration over 
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Fig. 15. CS J - 10-9 observations of L1448-C and the three 
NGC 1333 sources. In each plot the dashed line indicates the 
3cr detection limit. Lower 6 panels: Comparison between CS 
lines probing different excitation conditions in the envelopes, 
i.e., depths of NGC 1333-IRAS2 and -IRAS4B. The lowest ex¬ 
citation 3-2 lines are from IRAM 30 m observations the re¬ 
mainder from JCMT observations. 



Fig. 16. Comparison between CS lines probing different exci¬ 
tation conditions in the envelopes, i.e., depths of NGC 1333- 
IRAS2 and -IRAS4B. The lowest excitation 3-2 lines are from 
IRAM 30 m observations the remainder from JCMT observa¬ 
tions. 


+2 km s“^ from the systemic velocity: for the CS 7 = 10-9 
lines the emission integrated over this velocity range only con¬ 
tributes 20-25% of the total integrated emission and is still un¬ 
derestimated by the envelope models, especially for IRAS4B. 
Again this suggests that the observed CS J - 10-9 emission 


Table 11. Observed and predicted line intensities (f Tmb dv) 
for CS J — 10-9 for the three sources in NGC 1333. 


Source 

/mod [Kkm s^']“ 

Lbs [Kkms-'f 

IRAS2 

1.8 

< 2 

IRAS4A 

1.7 

12(3) 

IRAS4B 

1.1 

51(11) 


“Predicted CS J = 10-9 line intensity ado pting best fit abundances for 
each source from I.TOriJen sen e.t al.U2004cl) . *’Total line emission from 
Gaussian fits (IRAS4A and IRAS4B) or as 3(T upper limit (IRAS2). 
For IRAS4A and IRAS4B the number in parenthesis indicate the line 
emission integrated over ±2 km s“* from the systemic velocity. 


probes different material than traced in the bulk envelope ma¬ 
terial. 

Another simple estimate can be made assuming that the CS 
emission comes from a medium with a constant density and ki¬ 
netic te mperature. Using a non-LTE escape pro bability code, 
Radex jjansen et all Il994t ISchoier e~ ^ l2005l) . the CS col¬ 
umn density is estimated adopting a density of 3 x 1 0^ cm“^ 
and kinetic temperature of 100 K lIBlake et allll995li . consis¬ 
tent with the intensit y of the wing e i nission from the CS 5- 
4 and 7-6 lines from ljprgensen etd 1 (l2004cll . A high column 
density of ~ 5x10'^ cm“^ is needed to produce the observed CS 
10-9 emission. Even at such high column densities, the emis¬ 
sion is found to be optically thin. This column density is an 
order of magnitude larger th an found from the lower excitation 
lines bv iBlake et alJ (Il99-5li and corresponds to a; 5 - 10 % of 
the estimated CO abundance in the out flowing mate rial based 
on the CO column density calculated bv iBlake et al.L 

In contrast to the CS J — 10-9 lines, HDO 2n - 2 i 2 
is detected only toward NGC 1333-IRAS2 and not IRAS4A 
and IRAS4B (see Eig. I17> . Like CS 10-9 this line probes the 
warm, dense gas with an upper level energy of 90 K. The ob¬ 
served line is narrow (EWHM of a; 2.5 km s“') compared to 
the » 8 km s ' for the CS 7 = 10-9 and » 4 km s ' for the 
CH 3 OH lines toward IRAS4A and IRAS4B. This suggests that 
this line has its origin in a “hot inner region” of the NGC 1333- 
IRAS2 envelope, although relation t o the small-scale outflow 
(as seen in high resolution maps by Ijprgensen et al.t l2004bl) 
cannot be ruled out. Enhancements of HDO by up to a fac- 
tor of 10 were de rive d for the IRAS 162 93-2422 outflow by 
IStark et al.l (l2004l) . but IParise et°^ il200.5l) show that observa¬ 
tions of a larger number of HDO transitions are better explained 
by an abundance enhancement by two orders of magnitude in 
the innermost quiescent envelope. Observations of more tran¬ 
sitions will be needed before either scenario can be confirmed 
or ruled out for NGC 1333-IRAS2. Since both CH 3 OH and 
CH 3 CN observations indicate the presence of warm gas with 
abundance jumps in the inner (T < 90 K) region, NGC 1333- 
IRAS2 still seems the best candidate for further comparative 
studies of the passively heated, hot inner regions of low-mass 
protostellar envelopes. 
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Velocity [km s 

Fig. 17. Observations of the HDO 2ii - 2\2 line at 
241.5616 GHz toward the three NGC 1333 sources observed 
in the CH 3 OH 5^ setting. HDO is only detected toward 
NGC 1333-IRAS2. In each plot the dashed line indicates the 
3cr detection limit. 

4 . 3 . Comparison with IRAS 16293-2422 

The above results again raise the question what causes the 
richness of lines in IRAS 16293-2422 (even taking into its 
smaller distance and higher luminosity compared to, e.g., the 
NGC 1333 sources) since its abundances in the outer enve¬ 
lope are comparable to those of other sources. One explana¬ 
tion may be that its small-scale physical structure is signih- 
cantly different from that of the remaining sources. The central 
binary system has a separation of 8 - 10 " and affects the ma¬ 
terial in the envelope with emission from the various species 
centered around one or both components dMundv et all 1 199^ 
ISchoier et al.l l2004l) . Furthermore the circumbinary envelope 
appears to have an inner cavity of size comparable to the bi¬ 
nary separation. It may be that this relatively wide (~1000 AU) 
binary pushes material to larger distances where it can more 
easily be observed through single-dish observations with ~ 10 - 
15" beams. Alternatively, the circumbinary envelope may be 
heated on larger scales through each of the relatively luminous 
components compared with that expected from a single com¬ 
ponent in a simple spherical envelope. Finally, interferometer 
maps also show velocity gradients indicating that the outflow 
affects the envelope material close to the central protostar. It 
is possible that the outflow processing leads to abundance en¬ 
hancements of the organic species on larger scales than the pas¬ 
sively heated hot core, thus providing a larger hlling factor of 
the single-dish beam. 

5. Conclusions 

We have presented an analysis of H 2 CO and CH 3 OH line 
observations for a sample of 18 pre- and protostellar cores 
which have previously been studied through continuum and 
line observations and detaile d radiative transfer modeling 
( Ijdrgensen et all 12002 2004d). The se results complement the 
results by iMaret et alJ 11200411200.51) for a subset of sources 
In addition, observations and limits for high excitation CS 
7 = 10-9 transitions and lines of HDO and CH 3 CN are pre¬ 
sented for a subset of sources. Molecular abundances are de¬ 
rived through Monte Carlo line r adiative transfer and com pared 
to the results from the survey of LlOrgensen et alJ ll2004cli . The 
main conclusions are: 

- The H 2 CO data of most sources can be well-htted by con¬ 
stant abundances throughout their envelopes with an ortho¬ 


para ratio of 1.6+0.3. This implies thermalization of H 2 CO 
at low temperatures, e.g., on grain ice-mantles. Higher an¬ 
gular resolution data are needed to constrain the presence 
of any abundance jumps in the inner warm envelopes. 

- The H 2 CO abundances are related to the chemical network 
of the other species indicating that the same processes reg¬ 
ulate their abundances. As an example the H 2 CO data for 
NGC 1333-IRAS4A are well-ht by “drop abundance” pro- 
hles with a decrease in abundance from a few xl 0 “® - 10 “^ 
to a few xl 0 “^° in a limited part of the envelope, bounded 
inwards by the desorption temperature and outwards by a 
density corresponding to the lifetime of the core. A counter 
example is provided by NGC 1333-IRAS4B, where an 
abundance increase is only needed where the temperature 
rises above 20-30 K with no enhancement in the outermost 
regions. This indicates that for some sources other effects, 
such as the impact of an outflow, may be important for reg¬ 
ulating the H 2 CO abundances. 

- The upper limits to the CH 3 OH abundances for the entire 
sample are a fewxl0“'°-10“^. Th ese results are consistent 
with the abundances determined bv lRiickle & FulleTlll20o3i 
from lower excitation lines. 

- CH 3 OH observations for NGC 1333-IRAS2 and 
NGC 1333-IRAS4B require abundance jumps of at 
least two orders of magnitude at 90 K and 30 K, respec¬ 
tively. Together with the significantly broader lines of 
CH 3 OH compared with H 2 CO and other species, this 
suggests that the abundance increase for IRAS4B, in 
particular, is due to a compact outflow interacting with the 
nearby envelope. This is further supported by the broad 
high frequency CS 7 = 10-9 line, detected very strongly 
toward this source. HDO and CH 3 CN, in contrast, are only 
observed for NGC 1333-IRAS2 - possibly reflecting the 
presence of a passively heated, warm inner region for this 
source where molecules can evaporate. The CH 3 CN data 
for NGC 1333-IRAS2 also require a jump in abundance at 
90 K by about two orders of magnitude. 

This paper reinforces the importance of identifying (if pos¬ 
sible) unique chemical tracers of material heated “passively” 
by a central protostar and by shocked material in outflows. 
Even relatively high excitation lines from single-dish obser¬ 
vations (such as the CS 7 = 10-9 lines) may be affected 
by the outflow and can provide a good indication of the Ail¬ 
ing factor of dense shocked material. Future Herschel-HIFI ob¬ 
servations may provide additional tests of the chemical struc¬ 
ture by observations of high frequency lines, but the fact that 
shock and envelope chemistry are spatially unresolved may 
be problematic for these large beam data. Possibly the best 
way of distinguishing the different chemical scenarios will be 
through high angular resolution, high excitation observations 
with facilities such as the SMA and ALMA - or from stud¬ 
ies of outflows well- separated from the central protosta r as, 
for example, done iniBachiller & Perez Gutierrez! 1119971) and 
l.le)rgensen et al.l ll2004ah . Better knowledge about the physi¬ 
cal properties of the inner envelopes will also be important, 
since the adopted envelope models are extrapolations from the 
larger-scale observations of the cold dust in the outer enve- 
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lope. Infrared observations with the Spitzer Space Telescope 
can place better constraints there. 
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Appendix A: Line widths from gaussian fits 

Table IA.fl4A.2l lists the line widths (FWHM) for the observa¬ 
tions where Gaussians could be fitted. For the lower excitation 
H 2 CO observations not listed here the lines were integrated 
over +2 km s“*. 

Tab le A.l, Line widths (FWHM in km s ') for H 2 CO lines (see 
also lMaret et al.l ll2004l) '). 


Source 

Av(5o5 - 4o4) 

Av(5i5 -4i4) 

L1448-12 

1.5 

0.96 

L1448-C 

1.1 

b 

N1333-12 

1.9" 

_ b 

N1333-14A 

3.5" 

4.1" 

N1333-14B 

3.1" 

4.0" 

LI 527 

1.0 

0.9 

VLA1623 

0.63 

_b 

L483 

0.93 

1.4 

L723 

1.6 

2.3 

L1157 

0.52 

1.1 

CB244 

1.4 

1.8 

L1551 

1.6 

1.3 

L1489 


2.0" 

TMRl 


1.9 


Notes: “Non-Gaussian line profile. Intensity in Table Q refers to 
line integrated over ±2 km s“'. Line width defined as Av = 
f T dv/(1.064 Tpeak) where the integral is over the total line profile 
and Tpeak is the temperature at the line peak. ^Insufficient spectral res¬ 
olution to estimate width. 


Table A.2. CH 3 OH line widths (FWHM in km s *) for sources 
in NGC 1333. 


Line 

Frequency 1RAS2 

1RAS4A 

1RAS4B 

Sat - 4/f band; E-type 

-tOE 

241.7002 3.5 

6.0" 

3.4" 

-IE 

.7672 3.0" 

6.2" 

3.4" 

-4E 

.8132 4.9 



-t4E 

.8296 2.1 



-t3E 

.8430 

- 

-t 

-3E 

.8523 3.9 



-tlE 

.8790 3.4 

6.0 

3.4 

±2E 

.9044 3.4 

6.9 

3.5 

Sat - 4a: band; A-type 

-tOA 

241.7914 2.6" 

5.7" 

3.2" 

±4A 

.8065 3.2 



±3A 

.8329 3.4 

6.0 

2.8 

-2A 

.8430 

- 

-t 

+2A 

.8877 3.2 

6.0 

4.3 

7a - 6a band; E-type 

-IE 

338.3446 4.5" 

7.7" 

4.4'' 

-4E 

.5040 3.3 



-t4E 

.5302 3.2 



-3E 

.5599 3.0 



-h3E 

.5831 3.1 

5.3 

3.3 

-tlE 

.6150 3.3" 

7.0 

4.4 

7a - 6a band; A-type 

+0A 

338.4086 4.8" 

7.5" 

4.4" 

±4A/-2A 

.5127 3.3 

7.7 

5.4 

±3A 

.5419 4.6 

7.7 

5.6 

+2K 

.6399 3.0 

7.8 

4.3 


“t”The 5-4 -h 3E and -2A lines are blended at 241.8430 GHz. 
“Lines asymmetric. No Gaussian fitted. Line width defined as Av = 
/ Tdv/(1.064 Tpeak) where the integral is over the total line profile 
and Tpaak is the temperature at the line peak. 






















